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Abstract 

We propose a new scenario to produce the superheavy dark matter based 
on the inflationary universe. In our scenario, the inflaton couples to both a 
boson and a stable fermion. Although the fermion is produced by the inflaton 
decay after inflation, almost energy density of the inflaton is transmitted into 
the radiation by parametric resonance which causes the explosively copious 
production of the boson. We show that the fermion produced by the inflaton 
decay can be the superheavy dark matter, whose abundance in the present 
universe coincides with the critical density. We also present two explicit models 
as examples in which our scenario can be realized. One is the softly broken 
supersymmetric theory. The other is the “singlet majoron model” with an 
assumed neutrino mass matrix. The latter example can simultaneously explain 
the neutrino oscillation data and the observed baryon asymmetry in the present 
universe through the leptogenesis scenario. 
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The existence of the dark matter in the present universe is the commonly accepted 
consequence from observations [TJ. In addition, most of the inflation models [^]. which can 
solve the flatness and horizon problems, naturally predict the density parameter = 1. On 
the other hand, the big-bang nucleosynthesis implies that the contribution of baryons to the 
matter density at the present universe is at most 10%. Therefore, the present universe is 
dominantly fulfilled by the dark matter; Odm > 0.9. 

In general, there are two possibilities for the type of the dark matter. One is that the 
dark matter is a thermal relic, and the other is that it is a non-thermal relic. Most of 
discussions have been performed for the first type. In this case, the present abundance of 
the dark matter can be estimated, and, especially, we can derive an upper bound on its 
mass, mr>M < 500 TeV by the unitarity argument [||]. According to this argument, if the 
dark matter is superheavy [], it should be a non-thermal relic not to over-close the present 
universe. However, in this case, we should make it clear what is the mechanism to produce 
the superheavy dark matter in order to discuss its present abundance. 

Recently, some production mechanisms of the non-thermal superheavy dark matter were 
proposed ||, by which the dark matter is produced through gravitational effect, broad para¬ 
metric resonance and so on. Furthermore, some candidates for the superheavy dark matter 
were also considered in the context of the string theory, M-theory [|3J and supersymmetric 
theory with discrete gauge symmetry [||]. 

In this letter, we propose another scenario based on the inflationary universe in order 
to produce the superheavy dark matter whose abundance in the present universe coincides 
with the critical density. In our scenario, the inflaton very weakly couples to both a boson 
and a stable fermion whose masses are much smaller than the inflaton mass. Although the 
fermion is produced by the inflaton decay after inflation, almost energy density of the inflaton 
is rapidly transmitted into the radiation of the boson by parametric resonance [|7i] [|5| which 
causes the explosively copious production of the boson. If the boson couples to ordinary 
particles in the standard model with not so weak coupling constants, the universe can be 
thermalized as soon as the parametric resonance occurs. The production of the fermion is 
effectively over when the parametric resonance occurs. We will show that the fermion can 
be the superheavy dark matter, whose abundance in the present universe coincides with the 
critical value. 

We also present two explicit models as examples in which our scenario can be realized. 
One example is the softly broken supersymmetric theory. Note that it is natural for a boson 
to simultaneously couple to some bosons and fermions in the context of supersymmetric 
theories. The other is the “singlet majoron model” || with an assumed neutrino mass 
matrix. This example have other phenomenological and cosmological implications. In this 
model, we can simultaneously explain the solar and the atmospheric neutrino deficits and 
the baryon asymmetry in the present universe through the leptogenesis scenario 
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As mentioned above, in our scenario, the inflaton held couples to both the fermion and 
boson. Let us consider the inflaton potential of the form 


V=\\(f-a^ 


( 1 ) 


1 The word “superheavy” means that mass of the dark matter is larger than 500 TeV. 
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where 0 is the inflaton held, and it is regarded as a real held, for simplicity. We also consider 
interaction terms of the form, 

Ant = -9b4> 2 X ] X - , (2) 


where x and 0 are the boson and fermion helds which couple to the inhaton with (positive 
and real) coupling constants g B and gp, respectively. Assuming a -C m p i , where m p i ~ 10 19 
GeV is the Planck mass, the chaotic inflation occurs with the initial value of inhaton helds, 
(0) ~ (several) x m p i. Note that A ~ 10 -12 is implied by the anisotropy of the cosmic 
microwave background radiation |TT|| , and gp, Qf < vA is required by naturalness. 

Since vacuum lies at (0) = cr, let us shift the inhaton helds such as 0 —> 0 + a. Then, 
we can rewrite eqs.(|l]) and (^|) as 


A 


A , 3 A 

2 - ^ + 2 a0 + 8 < 


V = —a 2j ~ 2 ' ' - 4 


(3) 


and 


Ant = -^lO 2 + 2a0 + a 2 )x f X ~ 9f (0 + o-)00 , 


(4) 


respectively. Masses of the inhaton, boson x and fermion 0 are given by = y/\cr, 
rn 2 = g B cr 2 + m ( 2 and = g B a, respectively. Here, m 0 is tree level mass of the boson held. 
In the following discussion, we assume m x ,m^ -C m 

Now, let us consider behavior of an amplitude of the inhaton held after the end of 
inflation. In this epoch, the inhaton held is coherently oscillating, but its amplitude is 
decreasing due to both the expansion of the universe and the inhaton decay. When the 
amplitude becomes (0) -C a, the inhaton potential of eq. ([3]) is reduced to only the mass term, 
and thus we can regard the coherent oscillation as the harmonic oscillator with frequency 
0 = $ 0 cos(m^). 

Regarding the inhaton as the background held, equation of motion for the Fourier modes 
of x fi e ld is approximately given by 


dt 2 


Xk + (k 2 + 2g 2 B a^ 0 cos Xk = 0 , 


(5) 


where k is the spatial momentum, and we neglected mass term m 2 and g B (j) 2 term because 
of our assumption m x -C and (0) -C <7, respectively. Here, note that we also neglected a 
friction term 3 Hdxk/dt introduced by the expansion of the universe, and this approximation 
will be justihed later. The above equation is well known as the Mathieu equation [[1^] , and 
(narrow) parametric resonance occurs if the condition 4g B a^o -C m 2 is satished. In the 
following discussion, we take $ 0 ~ a and g B ~ A as a rough estimation. The amplitude of 
Xk with momentum in the hrst resonance band k ~ m^/ 2, whose contribution is dominant, 
is exponentially growing up such as Xk oc exp (/it), where 
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m<p 


~ m ^ . 


( 6 ) 


Note that the Hubble parameter is given by 
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H 


r^j 


(7) 


A a 2 $l 
m 2 pl 


1/2 





and H -C /i is satisfied because of the assumption a -C m p i. Then, the friction term 
3 Hdxk/dt can be neglected in eq.(||). Almost energy density of the inflaton field is rapidly 
transmitted into radiation of y fields by this parametric resonance. If y field (not so weakly) 
couples to ordinary fields in the standard model, the universe can be thermalized as soon as 
the parametric resonance occurs. Then, we expect that the reheating temperature ( T RP ) is 
estimated as 


V (</> 


cr 


' RH ■ 


( 8 ) 


Next, let us consider the energy density of the fermion which is produced from the 
time of the end of inflation {tEOi) till the time when the parametric resonance occurs and 
the universe is thermalized (t PR ). At t PR , the energy density of the produced fermion is 
estimated as 

P*{t PR ) = V(t EOI ) {l - e - r ^-^)} ’ ( 9 ) 

where V(tEoi ) ~ A m p i is the (potential) energy density of the inflaton at t E oi, a(t) is the 
scale factor of the universe, and T is the width of the inflaton decay process </> —> ipnp which 
is given by 


T = 



( 10 ) 


We regard time dependence of the scale factor between t E oi and t PR as same as the matter 
dominated universe 0, namely, a(t ) oc t 2 / 3 and H(t ) ~ (2/3)f _1 . The Hubble parameters, 
H{t E oi) and H(t PR ), are given by 


H(t E oi ) ~ ( V ^ ~ V^rnpi. , ( 11 ) 

V m P i ) 

( 12 ) 

V m h J m Pi 

If T H{t PR ) -C H (t eoi) 05 the energy density is approximately given by 

H{tpR) ~ (ffe/f) 

~ ^T H(t PR ) rn'h . (13) 

The present energy density of the fermion ip is given by 
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We can see that the condition is satisfied in our final result with our assumption <r <C m p i- 
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Pi>{t o) ~ Pip (t PR.) x 


T o_ 


3 


(14) 


Trh J 


x 


J 


where i 0 is the present time, and T 0 ~ 2 x 1CT 13 GeV is the present temperature of the 
microwave background radiation. If this energy density is comparable to the critical density 
(Pent ~ 8 x 10 _47 GeV 4 ), the fermion can be the dark matter. Using eqs. (I)-®, we can get 
resultant mass ratio of the dark matter to the inflaton Q. 

~ 10" 3 . (15) 


If we take, for example, a ~ 10 18 GeV, the fermion can be the superheavy dark matter with 
mass m^j ~ 10 9 GeV. 

Note that the existence of the parametric resonance is crucial in our scenario. If there is 
no parametric resonance, according to the old scenario of reheating, the reheating temper¬ 
ature is estimated as 

r = H ~ , (16) 

m p i 

where T is the decay width of the inflaton. However, in this scenario, the energy density of 
the fermion in the early universe is comparable with that of radiation, and thus the present 
universe is over-closed, if not gp is unnaturally small compared with gp. 

Here, we give a comment on our analysis in this letter. In general, there is a possibility 
that the broad resonance occurs before the amplitude of the inflaton is dumped into the 
narrow resonance band. However, it is non-trivial problem to definitely decide whether the 
broad resonance can affect or not, because of its stochastic behavior as analyzed in ||. Then, 
in our analysis we ignore the effect of the broad resonance for simplicity. If the energy of the 
inflaton is transmitted enough by the broad resonance, the period of reheating is shorten 
and total number of i/j produced by the inflaton decay is smaller than that of our estimation. 
In this case, the resultant mass of the superheavy dark matter becomes larger than that of 
our estimation in order for its present abundance to coincide with the critical value. It is 
always possible to understand our result as the lower bound on the super heavy dark matter 
mass. Some numerical calculations are needed for correct evaluation. 

In the following, we discuss two explicit models as examples in which our scenario of the 
superheavy dark matter production can be realized. The first example is the softly broken 
supersymmetric theory. Let us consider superpotential of the form 

W = X -m z Z 2 + X -m x X 2 + A z ZX 2 , (17) 

where Z and X are superfields with masses m z and m x , respectively, and X z is the 
dimension-less coupling constant. In the following discussion, we assume m x -C m z and 


3 We get the small coupling constant gp rs_/ 10 9 from eq.([l5|). It is easy to check that this coupling 
constant is too small for the superheavy dark matter to be in thermal equilibrium with other fields. 
This fact is consistent with our assumption. 
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Az -C 1, and identify the scalar component of Z and fermion component of X (i/ix) with 
the inflaton and the dark matter [] r respectively. The scalar potential is given by 

V = m 2 z Z f Z + m 2 x X f X + X 2 z (X f X) 2 + 4A|( Z f Z)(X f X) 

+ {A z m z ZX f2 + 2A z m x Z(X f X) + h.c.} . (18) 


Note that m z 
radiation [JTT[ . 


10 13 GeV is required by the anisotropy of the microwave background 


The inflaton held Z [] is coherently oscillating after the end of inflation and this oscillation 
is just the harmonic oscillator with frequency mz] Z = $0 sin (m z t). In the following, 
we consider the case in which the fourth term, 4A z (Z^Z)(X^X), in eq. (p~8D is dominant 
compared with the second line. For example, this case is realized, if we introduce the soft 
supersymmetry breaking terms which cancel the terms in the second line. Regarding the 
inflaton as the background held, the equation of motion for the Fourier mode of the scalar 
X is given by 


X" + ( A k - 2 q cos 2 z) X k = 0 , 


(19) 


where A k = k 2 /m 2 z + 2 q, q = (\ z $o/m z ) 2 , z = mt , and the prime denotes differential 
with respect to z. This is nothing but the Mathicu equation, and the (narrow) parametric 
resonance occurs if we hx parameters as \ z <C 1 and q ■ Cl, which means A^ <§C 10 -6 [] (we 
take $0 ~ m p i , see the following discussion). 

Since the oscillating inflaton is regarded as the harmonic oscillator from the beginning, 
the resonance occurs just after the end of inflation. Thus, it is naturally expected that 
H{tEoi ) ~ H(tpR ) and the reheating temperature is given by Tp„ ~ 

Following the discussion from eq.( 
component of X is described as 


z^o 




2 

P i- 


to eq. dT4|) , the present energy density of the fermion 




4 rn 2 / m z x 


TYlijjx 

Trh 


X 


T n 




Pcrit > 


( 20 ) 


where T = X 2 z m z /87T is the inflaton decay width with respect to the channel Z —> i^x^x- 
We can get the mass of the dark matter as m^ x 10 7 /A| 10 5 GeV. This result depends 

on the parameter X z . 

In the above discussion we assumed that the universe can be thermalized as soon as the 
parametric resonance occurs. This can be possible if we introduce the soft supersymmetry 
breaking term (A-term) such as 


C soft = AXHH , 


( 21 ) 


1 The fermion can be stable, if we introduce a discrete symmetry such as R-parity. 

5 We use the same notation both superfields and their scalar components in this letter. 

6 This coupling constant is also too small for the superheavy dark matter to be in thermal 
equilibrium with other fields. 
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where A is a parameter of mass dimension one, and H and H is the down- and up-type Higgs 
doublets in the supersymmetric standard model, respectively. Since the scalar X couples 
to the standard model particles by this soft terms, the thermalization of the universe is 
realized. 

Next, let us discuss the second example. There are some current data suggesting masses 
and flavor mixing of neutrinos. The solar neutrino deficit |H| and the atmospheric neutrino 
anomaly |T^| seem to be indirect evidences for the neutrino mass and flavor mixing from 
the viewpoint of neutrino oscillation. The “singlet majoron model” [§] is a simple extension 
of the standard model to give Majorana masses to neutrinos. It is well known that this 
model includes the see-saw mechanism [[HJ, which can naturally explain the smallness of the 
neutrino masses compared with other leptons and quarks. 

In addition, the model has a cosmological implication. Since, in the model, the lepton 
number is broken and CP is also violated by non-zero CP-phases in general, we can explain 
the observed baryon asymmetry in the present universe through the leptogenesis scenario 

0 - 


Assuming a typical matrix for the neutrino mass matrix, the singlet majoron model can 
be an example to which our scenario can apply. We can also simultaneously explain the 
observations of the solar and atmospheric neutrino deficits and the baryon asymmetry in 
the universe through the leptogenesis scenario. 

We introduce the Yukawa interaction of the form 


= -gyVL l 4>dV 3 R - g™v R ck (t)sV l R + h.c. , (22) 

where (j>d is the neutral component of the Higgs doublet in the standard model, (J) s is the 
electroweak singlet Higgs held, and i,j,k,l = 1,2,3 are generation indices. The Dirac and 
Majorana mass terms appear by non-zero vacuum expectation values of these Higgs holds. 
The neutrino mass matrix is given by 


0 m D 
m T M 

D J 


(23) 


where m D = g l -’(4>d) is the Dirac mass term, and M = g kl (4> s ) is the Majorana mass term. 
We assume typical Yukawa coupling constants and mass matrices as follows: 


'0 

a 

0 


'Mi 

0 

0 

0 

b 

0 

, M — 

0 

m 2 

0 

0 

b 

be lS 


0 

0 

M 3 


(24) 


Here, a, b and 5 are real parameters, and Mi, M 2 and M 3 are real masses of the heavy 
Majorana neutrinos. Note that, since the right-handed neutrino of the hrst generation 
decouples to the left-handed neutrinos, it can be stable if its mass is smaller than the singlet 
Higgs mass. This is the case we consider. 

In this model, we can identify the inhaton and the boson y in eqs.(jl|) and (|2) with the 
singlet Higgs (j) s and the standard model Higgs doublet, respectively. Then, our scenario can 
be realized, and thus the right-handed neutrino of the hrst generation can be the superheavy 
dark matter. The mass of the dark matter is three orders of magnitude smaller than the 
singlet Higgs mass, Mi ~ 10~ 3 x m^, according to our result of eq.(|l5|). 
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Next, we consider the problems of the solar and the atmospheric neutrino deficits. By 
the see-saw mechanism, the mass matrix of eq.(23) is approximately diagonalized as 


m n M 1 m T 0 

D D 

0 M ’ 


(25) 


where the mass matrix for the light neutrinos, m D M 1 m T D , is given by 


m„M 1 m T 

D D 


a 2 /M 2 ab/M 2 ab/M 2 

ab/M 2 b 2 /M 2 b 2 /M 2 

ab/M 2 b 2 /M 2 b 2 /M 2 + b 2 e 2iS /M 3 


(26) 


Note that Mi is absent in this matrix since the right-handed neutrino is decoupled to left- 
handed neutrinos. For simplicity, we assume a/b = e < Cl, M 2 /2M 3 e 7 < 1 and 5 -C 1. 
Neglecting the elements higher than the second order with respect to e, 7 and 5 , we can get 


b 2 


m n M m ~ — 
D D Mo 


e e 
1 1 
1 I + 27 


(27) 


The unitary matrix which diagonalizes this matrix is the physical Kobayashi-Maskawa 
matrix in the lepton sector at low energies. Moreover, assuming e C 7 , we can get hierar¬ 
chical mass eigenvalues such as b 2 /M 2 ( 0, 7 , 2 + 7 ) by diagonalizing this matrix. The mass 
squared differences and mixing angles corresponding to the solar and atmospheric neutrino 
oscillation data are given by 


A < ~ (ik) 72 '• s> " 2 2e@ ~ 4 ( 7 ) <<: 1 ’ 

f b 2 \ 2 

Ami ~ 4 \Tf 2 ) ; sin2 29 ® ~ 1 ' 


(28) 


By the assumption e < 7 , our description for the solar neutrino data corresponds to the 
small angle MSW solution |16[ . By appropriately choosing the values of the free parameters, 
we can reproduce the solar neutrino data Am 2 , ~ 5 x 10 6 eV 2 for the small angle MSW 
solution [|TTH and the atmospheric neutrino data A ml ~5 x 10^ 3 eV 2 




On the other hand, we can also explain the baryon asymmetry in the present universe 
through the leptogenesis scenario. Since the right-handed neutrinos in the second and the 
third generations couple to the left-handed neutrinos and m D of eq.([24|) includes CP-phase 


5, the original scenario of the leptogenesis can work. According to the original work |10 
can calculate the net lepton number production rate, which is proportional to 


we 


A = 


7TV 2 (7TlUn D ) 22 


Im 


mlm D )23('mUn D )l 2 


/I) 


(29) 


Here, /(M|/M|) ~ M 2 /2M 3 = 7 for M 2 -C M 3 , and v = 246 GeV is the vacuum expectation 
value of the standard model Higgs doublet. Using the Dirac mass matrix of eq.(24) and 


eq.(|28|), A is given by 
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In the following, to consider the observed baryon to photon ratio ns/n 7 ~ 10 10 , we use 
the rough estimation ns/n 7 r-*j A /g. 0, where g * ~ 100 is the effective degrees of freedom 
of the radiation in the early universe. 

Choosing appropriate values for the parameters M 2 and <5, we can simultaneously explain 
the correct abundance of the superheavy dark matter, neutrino oscillation data and the 
observed baryon ratio. When we take 6 = 0.1, for example, we can get the mass spectrum 
as follows: 

Mi ~ 10 9 < M 2 ~ 10 10 < M% ~ 10 11 < m<p 3 ~ 10 12 GeV . (31) 

Here, we took the vacuum expectation value of the singlet Higgs as a = 10 18 GeV. 

In summary, we proposed a new scenario to produce the superheavy dark matter. In our 
scenario, the inflaton couples to both the boson and the stable fermion whose masses are 
much smaller than the inflaton mass. The fermion is produced by the decay of inflaton after 
the end of inflation, but this production is effectively over when the parametric resonance 
occurs and the copious bosons are explosively produced. We showed that the abundance 
of the fermion in the present universe can be comparable to the critical value, and the 
fermion can be the superheavy dark matter. In addition, we presented two explicit models 
as examples in which our scenario can be realized. One is the softly broken supersymmetric 
theory. In this example, the dark matter has mass of larger than 10 5 GeV. The other 
example is the “singlet majoron model” with the assumed neutrino mass matrix. In this 
model, the right-handed neutrino of the first generation can be the superheavy dark matter 
with mass, for example, 10 9 GeV. Furthermore, this model can simultaneously explain the 
neutrino oscillation data and the observed baryon asymmetry in the universe through the 
leptogenesis scenario. 
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